Thermoelectric efficiency at maximum power in a quantum dot 
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We identify the operational conditions for maximum power of a nanothermoelectric engine con- 
sisting of a single quantum level embedded between two leads at different temperatures and chemical 
potentials. The corresponding thermodynamic efficiency agrees with the Curzon-Ahlborn expression 
up to quadratic terms in the gradients, supporting the thesis of universality beyond linear response. 
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The purpose of this letter is to present a detailed ther- 
modynamic analysis of electron transport through a sin- 
gle quantum dot connecting two leads at different tem- 
peratures and chemical potentials. Of particular interest 
to us is the efficiency of the thermal motor function, in 
which electrons are pumped upward in chemical poten- 
tial under the impetus of a downward temperature gra- 
dient. The study of this model addresses several issues of 
timely interest: nanotechnology, the study of thermody- 
namic properties of small devices that are prone to fluc- 
tuations, the question of universality for thermodynamic 
properties away from equilibrium, the role of quantum 
features in this respect, and the promise of thermoelec- 
tricity generated by nano-devices. We briefly comment 
on each of these topics. 

Spectacular technical and experimental progress in 
nano- and biotechnology have greatly increased our abil- 
ity to observe, manipulate, control, and even manufac- 
ture systems on a very small scale 0, |. In parallel, 
new theoretical tools and concepts have been developed 
that make it possible to exhibit the deeper relationship 
between fluctuations, entropy production and work, and 
the role of stochasticity in small scale non-linear nonequi- 
librium phenomena. In particular, stochastic thermody- 
namics formulates the thermodynamics of small entities 
subject to thermal fluctuations [E S Si Si 0- These 
developments are closely related to the celebrated fluctu- 
ation [8( and work theorems [9(. 

The concept of Carnot efficiency is a central corner- 
stone of thermodynamics. According to this principle, 
the efficiency, defined as the ratio of work output over 
heat input for a machine operating between two thermal 
baths at temperatures Ti and T r (T r > Ti) is at most 
equal to 
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case of Brownian motors, see [lfjj and references therein. 
However, reversible processes require infinitely slow op- 
eration, implying that such engines produce zero power. 
One of the important questions, when operating away 
from equilibrium, is the efficiency at maximum power. 
In a groundbreaking paper, Curzon and Ahlborn 
calculated this efficiency for the Carnot engine in the 
so-called endo-reversible approximation (taking into ac- 
count the dissipation only in the heat transfer process). 
They found a strikingly simple formula, namely 
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The equality is only reached for reversible operation. 
This is a universal result which remains valid for small 
scale fluctuating systems such as the well-documented 



Recently, it has been shown that the Curzon-Ahlborn effi- 
ciency is an exact consequence of linear irreversible ther- 
modynamics when operating under conditions of strong 
coupling between the heat flux and the work [hi . The 
value of 1/2 for the linear coefficient in Eq. ([2|) is there- 
fore universal for such systems. 

The efficiency at maximum power was also addressed 
in the context of stochastic thermodynamics in [l3j |. 
where it was shown that the efficiency at maximum 
power for a Brownian particle undergoing a Carnot cycle 
through the modulation of a harmonic potential is given 
by 77s = 2r? c /(4-r?c) « 7? c /2 + ^/8 + 3^/96 + .. .. By an 
entirely different calculation, dealing with the Feynman 
ratchet and pawl model (which operates under steady 
rather than cyclic conditions), the efficiency at maximum 
power was found to be [3| Vt = Vc/iVc — (1 — Vc) hi(l — 
Vc)] ~ Vc/2 + Vc/ 8 + 7?7c/ 96 + • ■ ■• All three of the above 
results agree, as they should to linear order in rj c . 
More surprisingly, the coefficient of the quadratic term is 
also identical. This raises the question as to whether uni- 
versality also applies to the coefficient of the quadratic 
term. 

The field of thermoelectricity went through a revival 
in the early 1990s due to the discovery of new ther- 
moelectric materials with significantly higher thermody- 
namic yields [HI]. Of particular interest are the devel- 
opments in the context of nanostructured materials [l|| . 
For example, thermoelectric experiments have been re- 
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FIG. 1: Sketch of the nanothermoelectric engine consisting 
of a single quantum level embedded between two leads at 
different temperatures and chemical potentials. We choose 
by convention T; < T r . Maximum power is observed in the 
regime e > fj,l > fir. 



ported on silicon nanowires (l7| . individual carbon nan- 
otubes (l9j and molecular junctions [l8| . Furthermore, it 
has been reported that Carnot efficiency can be reached 
for electron transport between two leads at different tem- 
peratures and chemical potentials, by connecting them 
through a channel sharply tuned at the energy for which 
the electron density is the same in both leads ^ 
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A double-barrier resonant tunneling structure has been 
proposed as a possible technological implementation [22] • 

The thermoelectric device whose properties we explore 
below is arguably the simplest prototype of such systems. 
It consists of a quantum dot with a single resonant energy 
level in contact with two thermal reservoirs at different 
temperatures, see Fig. Q] The dot can contain one single 
electron with a sharply defined energy e. The exchange 
of electrons between the leads through the dot will be 
described by a stochastic master equation [23|, 0, HH, 
and the corresponding thermodynamic properties can be 
obtained from stochastic thermodynamics 0, 0, 0, @, 0] • 
In anticipation of the forthcoming analysis, we note that 
this model displays perfect coupling between energy and 
matter flow: because of the sharply defined dot energy, 
every single electron carries exactly the same amount of 
energy. Hence, Carnot efficiency will be reached when 
operating close to equilibrium [lTJ, [2CJ, |21| , while Curzon- 
Ahlborn efficiency will be obtained at maximum power 
in the regime of linear response [13] • Going beyond these 
results, we will identify the operational conditions for 
working at maximum power. In particular, the efficiency 
at maximum power will be found to be r\ w t} c /2 + + 
. . ., with the coefficient of rfc again equal to 1/8. This 
provides further support for the thesis of universality for 
this value, especially since the regime of maximum power 
is found to lie entirely in the quantum regime. 

We now turn to the mathematical analysis of the ther- 
moelectric engine represented in Fig. [1] A single level 
quantum dot, with orbital energy e, exchanges electrons 
with a cold left lead, temperature Ti and chemical po- 
tential Hi, and with a hot right lead, temperature T r and 
chemical potential /i,.. The quantum dot is either empty 



(state 1) or filled (state 2). The crucial variables of the 
problem are the scaled energy barriers (with ks = 1) 



T ' 



v = Lr. 
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The exchange of electrons with the leads is described by 
the following quantum master equation 23|, 24, 25|: 
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The rates are given by 

W12 = £ w$ = J2 mi - u 

f—l, r v—l,r 

W21 = J2 w 2i = E a ^ 
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where /„ = [exp(x l /) + l] _1 is the Fermi distribution. The 
fact that a„ is independent of the dot energies is known 
as the wide band approximation. 

We are interested in the properties of the device at the 
steady state. The steady state distributions for the dot 
occupation follow from W 2 ipf s = W / i2P| s with pf s = 
1. The resulting probability current from the lead v to 
the dot is then 



(7) 



Using I r = —Xi and W12 + W21 = a r + ai,we can rewrite 
the result for the flux from the right lead as 



a(f r -fl), 



(8) 



where a = a r ai/(a r + ai). Eq. ([8]) is essentially the Lan- 
dauer formula for a single channel. 

The steady state heat per unit time Q r extracted from 
the lead r, and the steady state work per unit time 
(power) W performed by the device upon bringing elec- 
trons from right to left lead, are respectively given by: 

Qr = (e - Hr)lr = aT r x r (f r ~ fi) (9) 
W = (/i; - \i r )T T = aT r (x r - (1 - rj c )xi) (f r - /j) . 

(10) 

The corresponding thermodynamic efficiency reads 

? ?=7T = 7r = = i-( 1 - f ?c)— ■ (11) 

Qr Q r £ - Mr X r 

The entropy production associated with the master 
equation (g]) is given by [1, 0, % 0] 



<7 = £w£Vln 
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where i,j = 1,2. Noting that ln[W& ] /W%[ ] ] = 
finds, in agreement with standard irreversible thermo- 
dynamics [261 ] . the following expression for the entropy 
production: 

o = F m J m + F e J e = a(x t - x r )(f r - fi) > 0, (13) 

with thermodynamic forces for matter and energy flow, 
F m and F e , given by 



one Substitution of this result in (| 1 6|) gives 



F, 
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We stress that the corresponding matter and heat flow, 
given by 



Jm, — ^ri J p. — sF r 



(15) 



are proportional to each other. In other words, matter 
and heat flow are perfectly coupled and the condition for 
attaining both Carnot and Curzon-Ahlborn efficiency, 
namely, that the determinant of the corresponding 
Onsagcr matrix be zero, is fulfilled [3 El- 



We first discuss the case of equilibrium. Due to the 
perfect coupling, it is enough to stop one current, matter 
or energy, and the other one will automatically vanish. 
Under this condition, detailed balance is valid, X v = 0. 
It is clear from Eq. {5} that the matter flux (and hence 
also the energy flux) vanishes if and only if /; = f r or, 
equivalently, xi = x r . The efficiency then becomes equal 
to Carnot efficiency, cf. Eq. (jlip . and the entropy pro- 
duction vanishes, a = [cf. Eq. (|13|)]. Note that xi — x r 
does not require that the thermodynamic forces F m and 
F e vanish separately, i.e., at this singular balancing point 
equilibrium does not require temperature and chemical 
potential to be identical in both reservoirs 

[H El El El- 



We next turn to the operational condition for maxi- 
mum power. For given temperatures Tj and T r , we search 
for the values of the scaled electron energy barriers Xi and 
x r that maximize VV. From d Xl W = d Xr \V = 0, we find 
the following two equations determining these values: 



(fl - fr) + [Xr - (1 - Vc)xi] f r 



>k 



= o 

= 0. 



(16) 
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A first observation is that these equations depend only 
on the ratio of the two temperatures. Second, while the 
equations involve transcendental relations, one obtains 
the following explicit result by subtracting the first equa- 
tion from the second, 



xi = 2 In 



cosh{a; r /2} /cosh 2 {a; r /2} 



V 1 ~ Vc 



1 - Vc- 
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x r — \I~2 cosh(a; r /2) yj 2rj c — 1 + cosh(x r ) + 2(i] c — 1) 
cosh(x r /2) + \/2rj c — 1 + cosh(x^)/v2 
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cosh(a: r .) = 0. 



(19) 



Since an analytic solution of this equation is not possi- 
ble, we first turn to perturbative solutions for ry c close to 
the limiting values (reservoirs of equal temperatures) 
and 1 (cold reservoir at zero temperature). For the case 
rj c — > 0, we substitute x r — ao + a\T] c + Q2% + 0(Vc) 
in Eq. (|19[) and expand the resulting equation in rj c . 
The coefficients oi, a%, etc., are found recursively by 
solving order by order in rj c . At order zero, we find 
an identity. At first order, we find the transcendental 
equation ao = 2coth(ao/2). The numerical solution is 
ao = 2.39936. At second order and third order in r) c , we 
find a\ — — ao/4 and 02 = sinh(ao)/ {6 [1 — cosh(ao)]}. 
Substitution of these results and (fl8|) in (|11| leads to the 
following expansion of the efficiency at maximum power 
in the regime of small rj c : 



r? r,l [7 + csch 2 (q /2)] 
+ 96^ Vc 



(20) 



The expansion features the expected coefficient 1/2 for 
the linear term, but also supports the thesis that the 
coefficient of ?y 2 has a universal value, namely, 1/8. 

We next turn to the analysis of the case r/ c — ► 1 (TJ — > 
0). As we will see by self-consistency, the value of x r 
converges to a finite limiting value. With this a posteriori 
insight, we can easily identify the leading behavior of xi 
from Eq. (|18p . namely, xi ~ — ln(l — r) c ). Substitution 
of this result in (fl6|) leads to the conclusion that x r —> b, 
where b is a solution of the transcendental equation 
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(21) 



with numerical solution b — 1.27846. The corresponding 
efficiency converges to 1, albeit rather slowly, cf. Eq. (|ll|l . 

To complete the picture, we show in Fig. [5] the numer- 
ical solutions for xi and x r as functions of rj c as well as 
the corresponding rescaled maximum power W/(aT r ). 
Note that x r is always of order unity, so that the regime 
of maximum power cannot be well described by either 
a high or a low temperature expansion. It can also be 
seen that the maximum power is a monotonically increas- 
ing function of r\ c . The corresponding efficiency is repro- 
duced in Fig. [3] It is very close to the Curzon-Ahlborn 
efficiency, with relative deviations largest for large T] c . 

In conclusion, nanosystems with perfectly coupled 
fluxes, such as the quantum dot described here, are of 
great interest. They can operate as steady state Carnot 
engines. They probably possess universal features up 
to quadratic terms in nonlinear response when working 
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FIG. 2: (Color online) Scaled electron energy barriers xi and 
x r at maximum power, cf. Eqs. (|18|l and (|19|l . as a function 
of the Carnot efficiency r\ c — 1 — T\jT r . The dotted line 
represents the corresponding (scaled) power 10 x \V/(aT r ), 
cf. Eq. (TO). 




FIG. 3: (Color online) Efficiency at maximum power in a 
single level quantum dot as a function of the Carnot efficiency 
77c = 1 — Ti/T r (full line), as compared to Curzon-Ahlborn 
efficiency (dotted line) and Carnot efficiency (dashed line). 

at maximum power. One can speculate that they offer, 
from a technological point of view, advantages over their 
macroscopic counterparts. The above analysis can be re- 
peated for the quantum dot operating as a refrigerator, 
corresponding to the regime x r > x\. Such an analysis 
reveals that maximum cooling power (maximum Q r ex- 
tracted from the cold lead) is attained for x r — > oo and 
xi — > b, where b is again the solution to the transcen- 
dental equation (|21[) . The corresponding cooling power 
is thus bTi per transported particle, to be compared with 
the cooling power of 1) for a classical engine. 
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